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ABSTRACT: Mercury has perhaps the strangest behavior of any of the metals. Although 
the other metals in column IIB have an hep ground state, mercury's ground state is the body 
centered tetragonal /3Hg phase. The most common phase of mercury is the rhombohedral 
ctHg phase, which is stable from 79K to the melting point and meta-stable below 79K. 
Another rhombohedral phase, 7Hg, is believed to exist at low temperatures. First-principles 
calculations are used to study the energetics of the various phases of mercury. Even when 
partial spin-orbit effects are included, the calculations indicate that the hexagonal close 
packed structure is the ground state. It is suggested that a better treatment of the spin- 
orbit interaction might alter this result. 

INTRODUCTION 

Of all the metals in the periodic table, mercury has the most interesting behavior. A 
liquid at room temperature, the metal only solidifies below 234K. This phase, denoted aHg, 
(Pearson symbol hRl, space group R6m, Strukturbericht designation ALO), is a rhombohe- 
dral structure with one atom per unit cell and the primitive vectors at an angle of 7O°44.6'[0. 
Although ctHg has been seen experimentally down to 5K, below 79K the ground state is the 
phase /3Hg (prototype aPa, Pearson symbol tI2, space group I4/mmm, Strukturbericht des- 
ignation A a ), a body-centered tetragonal (bet) phase with a c/a ratio of 0.70710. Nearer 
to room temperature, aHg can be transformed to /3Hg by applying pressure. A third, 
meta-stable form, 7Hg, has also been observed 0. Like ctHg it is rhombohedral with one 
atom per unit cell, but the primitive vectors are at an angle of about 50°. 

Since these phases of mercury all involve only one atom per unit cell, first-principles ab 
initio methods are relatively easy to use. An extensive literature search, however, found 
very few studies of mercury]^, ||, [7|], with only one |7j including relativistic effects, which 
are important for all of the late fifth-row elements, and that only for metallic clusters. It 
is useful, therefore, to perform a series of first-principles calculations for the various phases 
of mercury. This paper presents results using the full potential, Linearized Augmented 
Plane Wave (LAPW) method|| |], |10| using the Hedin-Lundqvist[l 1 ] parametrization of 



the Local Density Approximation (LDA)|T^| to Density Functional Theory (DFT) |]T2"|, 13 1. 
The calculations were initially performed in the scalar-relativistic approximation ]14[, which 
essentially ignores the spin-orbit interaction while maintaining the remaining relativistic 
contributions. Spin-orbit corrections were then included using the "second-variational" 



method 15 



THE STRUCTURES OF MERCURY 

The primary structures of mercury are the rhombohedral aHg phase and the body- 
centered tetragonal /3Hg phase. Each phase can be described by two parameters: the 
volume and a parameter describing the orientation of the primitive vectors. In afflg this 
parameter is the angle a between the primitive vectors. In /3Hg the parameter is the c/a 



ratio of the tetragonal unit cell. Special values of these parameters lead to higher symmetry 
unit cells. The primitive vectors of the rhombohedral aHg phase can be written in the form 

a 1 = a(l+x, x , x ) 

a 2 = a ( x , 1 + x , x ) , x = §(y(l + 2 cosa)/(l — cos a) — 1) (1) 
a 3 = a(x,x,l + x) , 

where a is the angle between the primitive vectors. There are several special values of this 
angle. At a = the vectors (HD are collinear, while at a = 2ir/3 they are coplanar. These 
unphysical situations bound the range of a. Several high symmetry lattices can also be 
obtained from (TI]). At a = 7r/3 we find the fee lattice, at a = n/2 the simple cubic (sc) 
lattice, and at cos a = —1/3 the bee lattice. Because of these symmetries, a plot of the 
energy E(V, a) at fixed volume V would show the energy diverging as a approached both 
zero and 27r/3, with extremal points at a = n/3, 7r/2, and cos _1 (— 1/3). Since the aHg 
phase has a ~ 70° degrees, this phase will appear between the fee and sc phases. The 7Hg 
phase, with a ~ 50°, has a smaller angle than the fee phase. Figure |l] shows several of the 
phases found in the rhombohedral system. 
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Figure 1: Several of the structures found in the rhombohedral systemjl]. The a' phase shown 
here has an angle a between 7r/2 and cos _1 (— 1/3). The 7Hg phase would be to the right 
of the fee phase in this plot. 

The primitive vectors of the bet /3Hg phase can be written in the form 

a-x — ( a , , ) 

a 2 = ( , a , ) (2) 

a 3 — ( \& i \& i \c ) • 

This lattice is identical to the fee lattice when c/a = \/2, and to the bee lattice when 
c/a = 1. It is interesting to note that, within experimental error, the /3Hg lattice has 
c/a = 1/V2. At this value of c/a each mercury atom has two nearest neighbors, located 
directly above and below the atom along the z axis, at ±(2a3 — ai — a 2 ), and eight next- 
nearest neighbors. The mercury atoms thus form chains running along the z direction. 
Figure shows several phases in the bet system. 
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Figure 2: Several of the structures found in the bet system^. The j3' phase shown here has 
a value of c/a between the fee and bcc phases. 



FIRST-PRINCIPLES CALCULATIONS 

The first set of calculations were performed using the LAPW method in the "scalar- 
relativistic" approximation outlined above. This method essentially ignores the spin-orbit 
interaction, keeping the remainder of the relativistic corrections. In all of the calculations 
the muffin-tin radius was set at Rmt — 2.2 atomic units, a size chosen to allow large strains 
of the form ([!]) or (fj) while keeping the muffin tins from touching. The momentum cutoff 
Kmax was chosen so that Rmt K max = 10.5, yielding typical secular-equation dimensions of 
about 300 x 300. Increasing the cutoff to 11.5 decreases the energy by about 0.3m Ry for 
all structures and volumes. The K-point meshes were chosen using a regular mesh evenly 
spaced along the primitive vectors. Meshes of 150-200 K-points in the irreducible Brillouin 
zone, depending on the structure, yield total energies accurate to about 0.3 mRy compared 
to larger K-point meshes. The energies computed here are thus accurate to about 0.5 mRy. 

Comparing the scalar-relativistic and spin-orbit energies obtained from a Liberman- 
based atomic code[pl| show that the spin-orbit interaction contributes 29.5 mRy to the 
total energy. This is a relatively large contribution, so the spin-orbit interaction is included 
by a variational method[0. The spin-orbit energy is essentially converged if the second- 
variational basis uses 30 LAPW eigenstates in the variational calculation. 

Energy-volume curves were calculated for mercury in the fee, bcc and sc cubic lattices, 
the hexagonal close packed (hep) lattice, and the lattices described by ([I]) and (0) above. 
For convenience the latter structures will be referred to as aHg and /3Hg, even when they 
are outside the range of parameters which properly describe these structures. While the 
cubic structures require only knowledge of the volume to determine the structural energy, 
the hep, aHg, and /5Hg structures, require a knowledge of the energy as a function of the 
other lattice parameter. For the hep and /3Hg phases this parameter is c/a, while for aHg it 
is the angle a described in ([[]). Calculations were performed in both the scalar-relativistic 
approximation and with the variational spin-orbit energy included. 

The computations for aHg are shown in Figure [3. At most volumes there are two 
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Figure 3: LAPW calculation of the total energy of aHg in the rhombohedral structure (sys- 
tem!]). Both scalar- relativistic and spin-orbit results are shown. The error bars represent 
the estimated uncertainties in the energies. The lines between the points are drawn as an 
aid to the eye. 

minima in this plot, the global minimum between the sc and fee structures, corresponding 
to the observed aHg phase, and a secondary minimum between the bee and sc structures, 
which will be denoted c/Hg. These minima coalesce to the sc structure at large volumes. 
There are no minima in the region of the 7Hg phase @. Since this phase was obtained by 
shearing aHg, it is possible that the phase is stabilized by a non-hydrostatic shear. 

Similar calculations for /3Hg are presented in Figure |j. Again, there are two minima. 
The first, corresponding to the observed /3Hg structure, is near c/a ~ 0.7. The second, 
denoted /3'Hg, has a rather large shear of c/a « 1.7. This is analogous to the large c/a ratio 
found in the hep structures of Zn and Cd, and to the c/a ratio found in the calculations 
for hep Hg. Note that without the spin-orbit interaction the /3'Hg phase is lower in energy 
than the /3Hg phase, but the energy difference is not significant. The spin-orbit interaction 
lowers the energy of the /?Hg phase so that it is favored over the /3'Hg phase. 

The energy-volume curves for the low energy structures of mercury are calculated by 
finding the minimum energy as a function of the strain lattice parameter at each volume 
for the aflg, /3Hg, and hep Hg structures. The equilibrium energies for each of these phases 
is shown in Table [I], and the full energy volume curves are shown in Figure [|. From 
the calculations we must conclude that the hep structure is the ground state of mercury, 
contrary to experiment. The spin-orbit interaction does not change the relative ordering 
of the phases, and differences in the value of the spin-orbit interaction (the A column in 
Table |) are numerically insignificant. In addition, we see that the total energy calculations 
cannot distinguish between the aHg and the ct'Hg phases, nor between the /3Hg and the 
/3'Hg phases. 

If we nevertheless restrict ourselves to looking at the experimentally observed phases, 
we find the structural properties are in good agreement with experiment (Table |J). The 
calculated volume of the afig phase is about 13% smaller than the experimental volume, 




Figure 4: LAPW calculation of the total energy of /3Hg in the bet structure (system^). Both 
scalar-relativistic and spin-orbit results are shown. The error bars represent the estimated 
uncertainties in the energies. The lines between the points are drawn as an aid to the eye. 



while the /3Hg phase is only 1% smaller than the experimental volume. These results are 
consistent with LDA calculations for similar structures. The calculations do overestimate 
the angle (a) for the aHg phase, but get the c/a ratio correctly in the /3Hg phase. 



DISCUSSION 

Given the many successes of the DFT, and LDA in particular, in determining the 
structural properties of crystals, it is somewhat disturbing that we cannot correctly predict 
the ordering of the low-lying energy states of mercury. There are, of course, possible 
improvements to the LDA|l7|], and these may provide part of the answer. A more obvious 
problem with the present calculation is the form of the spin-orbit calculation. The present 



method||T0|, jig] uses scalar-relativistic orbitals as a basis for the second diagonalization of 
the Hamiltonian including the spin-orbit interaction. This is a good basis set for most 
states, but it has serious difficulties in dealing with the p states, since the relativistic pi/2 



Table I: The equilibrium energies for several structures of mercury, shifted so that the 
energy of the scalar-relativistic hep phase is set to zero. The second column shows the 
scalar-relativistic energy, the third the energy including the spin-orbit interaction, and the 
fourth the difference between the two. The "primed" phases are defined in the text. The 
energy of atomic Hg is also shown. 



Phase 


Scalar 


Spin- Orbit 


A 


Phase 


Scalar 


Spin-Orbit 


A 


fee 


.00244 


-.01013 


.01257 


a 


.00092 


-.01230 


.01322 


bee 


.00219 


-.01010 


.01229 


a' 


.00131 


-.01237 


.01368 


sc 


.00139 


-.01200 


.01339 


P 


.00196 


-.01160 


.01356 


hep 


.00000 


-.01309 


.01309 


F 


.00193 


-.01096 


.01289 


atom 


.02947 


-.00466 


.02947 
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Figure 5: LAPW calculation of the total energy Hg as a function of volume for the aHg, 
/3Hg, and hep phases. Both scalar-relativistic and spin-orbit results are shown. The error 
bars represent the estimated uncertainties in the energies. The lines between the points 
represent a Birch fit to the points shown. 



state is non-zero at the origin, while the scalar-relativistic p state vanishes there. This leads 
to an underestimation of the spin-orbit interaction of about 50% in this case, as can be 
seen in Table |[ In the atomic calculation, where the it is calculated exactly, the spin-orbit 
interaction contributes 29.5 mRy to the total energy. For the bulk structures, however, 
where the spin-orbit interaction is only approximated, it contributes about 13 mRy/atom 
to the total energy. It seems likely that a better treatment of the spin-orbit interaction 



will increase this contribution. It has been suggested fllPf that inclusion of pi/ 2 -like local 
orbitals in the LAPW basis would improve the spin-orbit energy, but this suggestion has 
yet to be implemented. 
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Table II: Structural properties of the aHg and /3Hg phases obtained from the LAPW 
calculations described in the text, compared to experiment. 





«Hg 








/3Hg 








V (a.u. 3 ) 


a (a.u.) 


a (deg) 


V (a.u. 3 ) 


a (a.u.) c (a.u.) 


c/a 


Experiment 


179.7 


5.643 


70.743 


152.1 


7.549 


5.338 


0.7071 


Scalar-relativistic 


158.2 


5.408 


75.6 


150.8 


7.52 


5.34 


0.71 


Spin-Orbit 


155.0 


5.372 


80.2 


150.6 


7.52 


5.33 


0.71 
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